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We investigated the magnetic phase diagram of the first Pr-based heavy fermion 
superconductor PrOs4Sbi2 by means of high-resolution dc magnetization measurements 
in low temperatures down to 0.06 K. The temperature dependence of the magnetization 
M(T) at 0.1 kOe exhibits two distinct anomalies at T c i = 1.83 K and T C 2 — 1.65 K, 
in agreement with the specific heat measurements at zero field. Increasing magnetic 
field H, both T C \(H) and T C 2(H) move toward lower temperatures without showing a 
tendency of intersecting to each other. Above 10 kOe, the transition at T C 2(H) appears 
to merge into a line of the peak effect which is observed near the upper critical field H C 2 
in the isothermal M(H) curves, suggesting a common origin for these two phenomena. 
The presence of the field-induced ordered phase (called phase A here) is confirmed 
for three principal directions above 40 kOe, with the anisotropic A-phase transition 
temperature T A : t| 100] > Tj 111] > T^ 10] . The present results are discussed on the 
basis of crystalline-electrical-field level schemes with a non-magnetic ground state, with 
emphasis on a Ti singlet as the possible ground state of Pr 3+ in PrOs4Sbi2- 
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1. Introduction 

The filled skutterudite compound PrOs4Sbi2 is considered to be the first Pr-based heavy 
fermion (hf) superconductor. The presence of the heavy quasiparticles has been suggested by 
the large specific heat jump AC/~/T c = 500 mJ/K 2 mol at T c = 1.85 K, 1 ) and more directly by the 
recent de Haas-van Alphen experiments. 2 ) PrOs4Sbi2 has received considerable attention, not only 
because it is the first Pr-based hf superconductor but also due to its unusual superconducting (sc) 
properties. Absence of a coherence peak in the temperature variation of the nuclear spin-lattice re- 
laxation rate in the Sb-NQR experiment is suggestive of non s-wave superconductivity. 3 ) However, 
it is under debate whether the sc gap has nodes or not. The thermal conductivity measurement in 
rotating magnetic fields has reported that the sc energy gap has point nodes. 5 ) On the contrary, 
the muon-spin relaxation 4 ) and the NQR 3 ) experiments on PrOs 4 Sbi2 have reported an isotropic 
superconducting (sc) energy gap. 

Unconventional superconductivity in PrOs4Sbi 2 is also suggested by its unusual phase diagram. 
The recent specific heat measurements at zero field revealed double sc transitions at T& = 1.75 K 
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and T C 2 = 1.85 K. 6,7 ) Moreover, the thermal conductivity measurement has claimed that the gap 
symmetry changes from two-fold to four-fold as the field increases. 5 ) Both experiments suggest the 
presence of the multiple sc phase in PrC^Sb]^. 5-7 - 1 However, this interesting possibility certainly 
needs to be verified by various experimental methods. 

One of the important issues in PrOs4Sbi2 is how the hf state is formed. In this regard, a two- 
channel quadrupolar Kondo effect has been proposed on the assumption that the CEF ground state 
for trivalent Pr ion is a non-Kramers T3 doublet. 8 ) On the contrary, Aoki et al. reported that the 
specific heat data of PrOs4Sbi2 are better described by a Fi ground state model. 9 ) Thus, the CEF 
ground state in PrOs4Sbi2 has not been established yet. In this respect, it would be important to 
notice that the recent specific heat measurements on PrOs4Sbi2 revealed a field-induced ordered 
phase in high fields above 40 kOe for the [100] direction. 9 ) Since such a field-induced ordering has 
never been observed in the Ce-based hf systems, clarification of its origin would be a clue to resolve 
the electronic state in this compound. 

From the above point of view, we investigated the sc state and the magnetic phase diagram 
of PrOs4Sbi2 in detail by dc magnetization measurements. The temperature dependence of the 
magnetization M(T) is found to exhibit the double sc transition in magnetic fields up to 14 kOe. 
In particular, the presence of the A-phase is found not only for [100] direction but also for all the 
principal directions. 

The paper is organized as follows. After a brief description of the experiment in Sec. II, we 
show the magnetization results in Sec. III. In Sec. IV an analysis based on the CEF model is given, 
followed by a discussion on the sc phase and phase A in Sec. V. The paper closes with a summary 
in Sec. VI. 

2. Experimental 

A single crystal of PrOs4Sbi2 (20 mg) was grown by the Sb-self-flux method, 10 ' 11 ) using high- 
purity raw materials 4N (99.99% pure)-Pr, 4N-Os and 6N-Sb. For the dc magnetization measure- 
ment below 2 K, we used a capacitive Faraday force magnetometer installed in a 3 He- 4 He dilution 
refrigerator. 12 ) The resolution of the system is better than 10 -5 emu. Throughout the measure- 
ments, we applied a field gradient (1 kOe/cm) which is indispensable to this method. The resulting 
field distribution inside the sample was estimated to be less than 150 Oe. The magnetic field was 
applied along the three principal directions of the body centered cubic PrOs4Sbi2 structure: [100], 
[110], and [111]. A SQUID magnetometer (MPMS, Quantum Design Co.) was also used to measure 
the dc magnetization in the temperature range of 2-300 K in field up to 70 kOe. 
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Fig. 1. Magnetic susceptibility M(T)/H(= x( T )) of PrOs 4 Sbi 2 in a magnetic field of 2 kOe applied along the [100] 
direction. The arrow indicates the superconducting transition temperature T c \. The inset demonstrates the double 
transition at T c i = 1.51 K and T c2 = 1.31 K in the M(T)/H curve for H || [110] at H = 6 kOe. 

3. Results 

Figure 1 shows the temperature dependence of the magnetic susceptibility M(T)/H(= x(T)) of 
PrOs4Sbi2 in a fixed field of 2 kOe applied along the [100] direction, plotted on a semidogarithmic 
scale. Above 50 K, the x(T) curve follows a Curie- Weiss law with the effective moment /x e fr=3.5 /j,b 
and the Weiss temperature #cw=-15 K. The obtained values of /j, e g and #cw are in good agreement 
with the previously reported ones. 1 ) The x(T) curve shows a rather broad peak around 3.5 K as 
reported before, 1 ) though the origin is not clarified yet. On further cooling, \ continues decreasing 
until a sharp kink due to the superconductivity appears at 1.51 K. 
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Fig. 2. Temperature derivative of the magnetization dM/dT plotted as a function of T at several fields below 14 kOe 
for two field orientations. The data were taken upon gradually warming the sample after field cooling. Two types 
of arrow indicate the double transition at T c i and T C 2- 



Double sc transition at T c \ 



= 1.85 K and T C 2=1.75 K in zero field has been reported in the 
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specific heat C(T) measurements of PrC^Sb^. 6,7 -* To check the presence of the double transition, 
we carefully measured the temperature dependence of the magnetization in low fields below 14 kOe. 
Temperature derivatives of the magnetization dM(T)/dT at several fields for H | [100] and [110] 
are shown in Fig. 2. The data were taken upon gradually warming the sample after cooling in 
magnetic field from above T c \. In the dM(T)/dT curve at H = 100 Oe, we found double transition 
at T c i = 1.83 K and T C 2 = 1.65 K, which are close to those determined in the C(T) measurements 
at zero field. In the C(T) measurements the double transition was not visible at 5 kOe. 6 ' 7 ) Here 
we were able to trace the two anomalies in dM/dT up to 14 kOe, as seen in Fig. 2. With increasing 
the field, both anomalies shift toward the lower temperature side. Interestingly, the anomaly at 
T c i becomes weaker, whereas the one at T C 2 becomes more pronounced. The inset of Fig. 1 shows 
the double transition at 6 kOe in the M{T)/H curve for H || [110], in which a small but distinct 
magnetization jump can be seen at T C 2. 




H(10kOe) 

Fig. 3. Isothermal magnetization curves of PrOs4Sbi2 at the base temperature of 0.06 K for three principal directions 
[100], [110], and [111]. The inset shows the M(H) curve for [100] and [110] direction around the peak effect. Two 
distinct anomalies at H' and H C 2 are shown in the inset. 

The magnetization process M(H) at 1.9 K, which is well above the ordering temperatures, 
are shown in the upper inset of Fig. 3. The magnetization is nearly isotropic. The difference 
between the M values for two directions is only 5 % even at the high field of 125 kOe, with the easy 
magnetic direction being [110]. The M(H) curves of PrOs4Sbi2 at the base temperature of 60 mK 
for three principal directions [100], [110], and [111] are shown in Fig. 3. The magnetic anisotropy 
remains small even in the base temperature of 0.06 K. In the M(H) curves, superconductivity 
manifests itself in the hysteresis at low fields. A marked peak effect is seen around 15 kOe for all 
the directions. As seen in the lower inset of Fig. 3, the peak effect is considerably anisotropic; the 
peak effect in the [110] direction becomes most pronounced. 

In Fig. 4 we plot the differential susceptibility dM/dH for the increasing-field process for three 
principal directions. The dM/dH curve for [100] direction at 0.06 K in the vicinity of H C 2 is shown 
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in the inset of Fig. 4. Here we define the upper critical field H C 2 by a sharp kink in the dM/dH 
curve as indicated in Fig. 4. The so-derived H° c2 = H c2 (T = 0.06 K) for H || [100] is 22.4 kOe, 
which coincides well with H® 2 determined from the resistivity measurements. ^ H® 2 for H || [100] is 
1.4 % higher than that (=22.1 kOe) for H || [110] and [111]. When the field is increased, the sharp 
structures due to the peak effect appear below H&. The onset of the peak effect is characterized by 
a dip in dM/dH, as usual. Interestingly, the peak effect abruptly diminishes at the field denoted 
by H' where a discontinuous jump is observed in dM/dH. This is a new feature that has not been 
observed in the conventional peak effect, suggesting a certain dramatic change in the vortex state. 
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Fig. 4. Differential susceptibility dM(H)/dH for the increasing field process at several temperatures below 0.88 K. 
The results for the three principal directions [100], [110], and [111] are shown. For clarity, the origin of each curve 
is shifted vertically. The arrows indicate various transition fields. The inset shows the dM/dH in the vicinity of 
H c2 for H || [100]. 

Going back to the data in Fig. 3, a clear upward bend of M{H) is visible around 45 kOe 
in the normal state for all the directions, indicating that a second-order phase transition sets in. 
The present result is consistent with the previous C (T) measurements which revealed the presence 
of the field- induced ordered phase (called phase A here) for H || [100]. 9 -* The dM/dH curves in 
Fig. 4 show the A-phase transition more clearly. For H || [100] only a single anomaly is visible 
at i?A=43 kOe above H C 2, which is defined as the transition field from the normal state to the 
A-phase. This anomaly gradually shifts towards the high-field side as the temperature increases. 
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For H || [110] four distinct anomalies are found at 47, 81, 97, and 116 kOe in the dM/dH curve 
at 0.06 K. From the lower field side we assign the four critical fields as Ha, Hi, H2, and H' A , 
respectively. H' A is considered to be the field where the A-phase disappears, whereas H\ and H2 
probably come from certain change in the ordered structure. As the anomaly at H 2 is markedly 
sharp, this phase transition is likely of first order. With increasing the temperature, Hi and H2 do 
not shift very much, whereas Ha and H' A approach to each other gradually. For H || [111] there are 
two anomalies at i?A=45 kOe and H' k =Yl\ kOe in the dM/dH curve at 0.06 K. These anomalies 
likewise approach to each other as the temperature increases. The de Haas-van Alphen (dHvA) 
effect is clearly seen above 80 kOe for all the directions, which is indicative of the excellent quality 
of our sample. The dHvA frequency is estimated to be about 1.03 x 10 3 T, which coincides well 
with that of the (3 branch obtained by the field modulation technique. 2 ) 
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Fig. 5. Temperature dependence of the M/H value at several fields above 50 kOe for three principal directions [100], 
[110], and [111]. 




Fig. 6. Temperature derivative of the magnetization AM /AT at several fields above 50 kOe for three different direc- 
tions [100], [110], and [111]. The arrows indicate the A-phase transition temperature Ta. 
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We next move on to the temperature dependence of the magnetization M/H. Figure 5 shows 
the M{H)/H curves at several fields above 50 kOe for three principal directions. As already shown 
in Fig. 1, M/H at low field exhibits a broad maximum around 3.5 K. This feature can still be seen 
in the M/H data for 50 kOe which show a strong increase above 1 K. One can see the A-phase 
transition as a sudden change of the slope in the M{T)/H curve at around 0.7 K. To show the A- 
phase transition more clearly, the temperature derivatives of the magnetization dM/ dT are shown 
in Fig. 6. We can easily determine the A-phase transition temperature Ta by the position of a step 
in the dM /dT data as shown by the arrows. Ta first increases with H, taking the maximum at 
around 90 kOe, and then decreases at higher fields. Although not clearly seen in Fig. 5 due to a 
limited temperature window of the present measurements, the broad maximum in M/H appearing 
at low field (Fig. 1) gradually shifts to the low field side as H increases. The maximum eventually 
disappears at around H ~ 90 kOe, above which M/H becomes a decreasing function of T. It 
should be noted that the field where the broad peak in the M(T) curve disappears is close to the 
field where Ta is the highest. The observation suggests that the origin of the maximum in x is 
closely related to the existence of phase A. 
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Fig. 7. H — T phase diagram of PrOs4Sbi2. A, P, and SC represent for the A phase, paramagnetic phase, and the 
superconducting phase, respectively. Open and closed symbols were determined by the dM(T)/dT and dM(H)/dH 
data, respectively. Since the superconducting phase is nearly isotropic on this scale, only the data for H || [100] is 
shown in the figure. Cross symbols denote to the onset of the peak effect in the M(H) process. 

The present magnetization data are summarized in the H — T phase diagram displayed in 
Fig. 7. Open and closed symbols are the transition points derived from the M(T) and M{H) 
data, respectively. T ci and T C 2 are nearly isotropic. Since T C 2(H) smoothly connects to H', the 
origin of T C 2 is presumably the same as the peak effect observed in the M{H) curve. Recently, a 
multiple sc phase diagram has also been reported from the thermal conductivity measurements on 
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PrOs4Sbi2- 5 - ) They observed that a gap symmetry changes from two-fold to four-fold at a certain 
field. Its boundary, beginning near T C 2(0), increases up to 0.7 T at about 0.5 K with decreasing T. 
The boundary T C 2(H) determined in the present experiment is, however, essentially different from 
that reported in ref. 5 ) So far, we have not succeeded in finding any appreciable anomaly in M(H) 
around 0.7 T as reported in ref. 5 ) 

Regarding the field-induced ordering, it should be emphasized that the A-phase exists for all 
the directions. Similar results have also been reported independently in ref. 13 ) We found that the 
anisotropy of the ordering temperature Ta is T^ 00 \h) > T A m '(iT) > T^ 10 \h). Due to the unique 
phase diagram, it has been proposed that the A-phase is an antiferro-quadrupolar (AFQ) state. 7 ' 9 ) 
While this is feasible, however, we cannot identify the order parameter only from the shape of the 
phase diagram, as will be shown later. We note that the anisotropy of Ta in PrOs4Sbi2 is different 
from that of the AFQ transition temperature Tq in the O xy -type AFQ compound CeBe 14 ** and the 
O^-type AFQ compound PrPb 3 . 15 ) 

4. Analysis 

In this section we analyze the thermodynamic properties of PrOs4Sbi2 by a CEF model. We 
believe that the localized /-electron picture (4/ 2 ) is a good starting point to understand the mag- 
netic properties of the system, 2 ) though the hybridization effect between the / electrons and the 
conduction electrons is surely important. We introduce the following mean field Hamiltonian: 



H = H C EF-g^BJ-(H + n(M)), (1) 

where the gj is the Lande g , -factor(=4/5) and is the Bohr magneton. Hcef is the CEF Hamil- 
tonian for the subspace of the J = 4 multiplet in the Oh group, 16 ) and can be written as 

n° — ^n 4 n° — 91 d 4 

ffCEP ^ [I 5i_^ +(1 _ N) 2^i 1 , (2) 

where the parameter W characterizes the overall strength of the CEF potential at the Pr site and 
x parameterizes the relative strength of the sixth-order Stevens' operators, n is a coefficient for the 
mean field magnetic exchange interaction. 

Bauer et al. proposed two CEF energy level schemes to explain the observed peak in the 
magnetic susceptibility x(T) at 3.5 K. 1 ) Parameters of the CEF interaction for the two schemes 
(referred to as schemes A and B) are shown in Table I. In Table II the energy levels for these 
schemes are given. An important difference between these two level schemes is the ground state. 
The ground state in scheme A is a non-Kramers r% doublet, whereas that in scheme B a T\ singlet. 
The first excited state in both cases is the same magnetic Ts triplet, located at 11 K and 6 K for 
schemes A and B, respectively. Since the second and third excited levels in both schemes have 
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much higher energies, the low-temperature and low-field properties are dominated by the low-lying 
two levels. 

Table I. Parameters of the CEF interaction and the magnetic molecular field coefficient for the two different CEF 
level schemes. 



W (K) x n (g/emu) 



scheme A -5.44 -0.72 -2300 
scheme B 1.85 0.5 -6000 



Table If. Energy levels for the scheme A and B 



scheme A 


scheme B 


Ti(313 K) 


r 3 (in k) 


r 4 (130 K) 


T 4 (65 K) 


r 5 (n K) 


r 5 (6 k) 


r 3 (o k) 


Ti(0 K) 



Let us first discuss the present magnetic susceptibility data x(T). In a low field limit, the 4/ 
part of the susceptibility obeys 

— = — - n, (3) 

X4/ XCEF 

where xcef denotes the single- ion magnetic susceptibility. The calculated %4/(T) for schemes A 
and B are shown in Fig. 8. The coefficient n was determined so that x&f best reproduces the 
result at 2 kOe (solid line). The obtained values of n for two schemes are shown in Table I. A 
main difference between the two calculated x&f 1S 1X1 the low-temperature behavior below 4 K. The 
X4f(T) curve for scheme A becomes nearly constant just after showing a weak peak at 4 K. On 
the other hand, Xif{T) for the scheme B strongly reduces upon cooling below 4 K, because of a 
difference in the Van Vleck contribution between the ground state and first excited state. Clearly, 
the XAf curve for scheme B is in better agreement with the experimental result, which continues 
decreasing down to T c = 1.5 K. Of course, we cannot judge the ground state of PrOs4Sbi2 from 
this fact alone, since the hybridization effect cannot be neglected in this system and the simple 
localized electron model might not be applicable at such a low-temperature and low-field region. 
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Fig. 8. Magnetic susceptibility of PrOs4Sbi2 for H | [100]. The dashed and dotted lines are the calculated x(T) for 
scheme A and B, respectively. Experimental data at H — 2 kOe (solid line) are also shown for comparison. 




Fig. 9. Field dependence of the two lowest energy levels. Upper panels show the Zeeman effect for F3 ground state 
model (scheme A), whereas the lower ones are those for the Y\ ground state model (scheme B). 
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Fig. 10. Magnetization curves M(H) at the paramagnetic state calculated at T — 2 K for the T3 (scheme A) and 
the Ti (scheme B) ground state models. 

We next calculate the Zeeman splitting of the CEF levels for the two orientations [100] and 
[110] by using the parameters in Table I, and the results for the relevant low- lying levels are given 
in Fig. 9. In scheme A, the Zeeman effect is anisotropic. Accordingly, the calculated magnetization 
curves are strongly anisotropic, as shown in Fig. 10. It should be noticed that a level crossing in 
the ground state appears at 190 kOe only for H || [100], as shown in Fig. 9. On the contrary, the 
field dependence of the CEF levels in scheme B is almost isotropic, leading to a nearly isotropic 
magnetization curves as shown in Fig. 10. It should be emphasized that a level crossing occurs at 
60 kOe for both orientations, though the crossing for [110] is weakly of repulsive type. 

Aoki et al. measured the specific heat of PrOs4Sbi2 in magnetic fields and obtained the entropy 
S as a function of T for various fields. 9 ) They revealed that the entropy increases with H at low 
temperatures below ~ 3 K. The observation is rather unusual because the entropy is normally a 
decreasing function of magnetic field. They argued that the observed decrease in the entropy by 
the field is related to the CEF level crossing expected in scheme B. Here we extend the mean-field 
Hamiltonian (1) to a two-sublattice model to treat the field- induced ordering, and calculate the 
temperature dependence of the entropy S(T) at various fields below 60 kOe for schemes A and B 
(Fig. 11). The calculated S(T) for scheme A decreases in the whole temperature range as the field 
increases, whereas for scheme B the S(T) value below ~ 2 K increases with the field, which is in 
agreement with the experimental results by Aoki et al. 9 ^ The kink at T*=0.7 K in the S(T) curve 
for 60 kOe is due to an AF magnetic transition caused around the level crossing. It can be shown 
that the field-induced AF ordering appears for all the field directions because of the nearly isotropic 
Zeeman effect. An AFQ transition is equally possible if we include an AFQ interaction in eq. (1). 
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Fig. 11. Temperature dependence of the calculated entropy S(T) at several fields along the [001] direction for the 
Ta (scheme A) and the Fi (scheme B) ground state models. 

On the other hand, it is obvious that no AF ordering can be expected for the scheme A (upper 
panel of Fig. 11) with the given weak AF interaction, since the CEF ground state is non-magnetic 
and the level crossing occurs only for H || [100] at a very high field (~ 200 kOe). Even if we take a 
quadrupolar interaction into account, it is very hard to explain the field-induced ordered phase (cf. 
Fig. 7) with the T3 ground state model, so long as we stay within the localized electron picture. 
In this situation, an AFQ phase should appear from zero field because the CEF ground state itself 
possesses the quadrupolar degree of freedom. This AFQ phase may exhibit a reentrant behavior 
like the case in PrPb 3 , 15) but never exhibits a field-induced ordering like the one shown in Fig. 7. 

5. Discussion 

In the preceding section, we showed within the localized electron picture that the magnetic 
properties of PrOs4Sbi2 are better explained by the T\ ground state model. In particular, there is 
a serious difficulty in the T% ground state model on explaining the field-induced ordering. Of course, 
we cannot go too far with the localized electron picture because the 4/ electrons in this compound 
are considered to be delocalized at low temperatures due to hybridization with conduction electrons. 
If the delocalization effect is strong enough to destroy the quadrupolar ordering at low fields, then 
a completely different situation might appear in the T3 ground state model: a scenario that a 
quadrupolar ordering is restored in strong field by suppression of the delocalization effect. A 
magnetic analogue of this scenario might be a field-induced antiferromagnetic phase transition in 
a system where antiferromagnetic interactions among localized magnetic moments are competing 
with the Kondo effect. The fact that such a transition has not been observed yet in any Ce-based 



Magnetic Phase Diagram of the Heavy- Fermion Superconductor PrOs4Sbi2 13 

hf compounds suggests that the aforementioned scenario for T3 is unlikely. 

A remaining problem with respect to the Ti ground state model is how the heavy-fermion 
state is formed at low temperatures. Concerning this, there are several theoretical studies, which 
attempt to explain the Fermi liquid state with heavy quasiparticles in the f 2 configuration with a 
singlet CEF ground state by taking the CEF excited states into consideration. 18 ' 19 ) Actually, many 
U-based-hf compounds (e.g. UPt3, and URu2Si2), which are likely to be in the f 2 configuration, is 
believed to possess a singlet CEF ground state. Accordingly, formation of the hf state of PrOs4Sbi2 
might be understood by the Ti ground state model. 

As for the double sc transition, at present there is no clear-cut evidence that this phenomenon 
is intrinsic. Since the phase boundaries T c i{H) and T C 2(H) are likely to be scaled to each other 
(Fig. 7), there remains a possibility of a phase separation in spite of the high quality of the sample: a 
small fraction of the volume first undergoes a sc transition at T cl , and at slightly lower temperature 
T C 2 the whole volume of the sample changes into the sc state. To eliminate this scenario, further 
experimental efforts will be needed. 

6. Summary 

We measured the dc magnetization of PrOs4Sbi2 in low temperatures down to 60 mK. The 
presence of the double sc transition, which was first observed in the specific heat measurements at 
zero field, is confirmed in fields up to 14 kOe. The field dependence of T C 2 is qualitatively the same as 
that of T c i. We also observed phase A not only for [100] direction but also [110] and [111] directions. 
The anisotropy of the A-phase transition temperature is T^ 00 ^ > T^ 111 ' > T^ 110 '. For H || [110] two 
additional transitions are found inside phase A. Comparison between the experimental results and 
the calculations on the basis of the CEF model suggests that the CEF ground state of Pr 3+ is a T± 
singlet, in contrast to the previous report. 1 ' 7 ) 
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